Jet quenching at intermediate RHIC energies 
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The final state energy loss of fast partons penetrating a longitudinally expanding quark-gluon 
plasma of effective gluon rapidity density dN 9 /dy = 650 — 800 is evaluated and incorporated together 
with the multiple initial state Cronin scattering in the lowest order perturbative QCD hadron 
production formalism. Predictions for the neutral pion attenuation in central Au + Au collisions at 
the intermediate RHIC energy of ^/snn = 62 GeV relative to the binary collision scaled p + p result 
are given. The quenching is found to be a factor of 2 — 3 with a moderate transverse momentum 
dependence and the attenuation of the away side di-hadron correlation function is estimated to be 
3-5 fold. 
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I. INTRODUCTION 

Recent combined experimental measurements of the 
nuclear modification to the moderate and large trans- 
verse momentum hadron production in Au + Au Q, Q 
and d + Au [3] reactions have provided strong evi- 
dence in support of the dominance of multiple final 
state interactions Q over the initial state Cronin scat- 
tering 5] and possible nuclear wavefunction effects |6J 
in relativistic heavy ion collisions. These findings pave 
the way for detailed studies of derivative jet quench- 
ing observables 0, @ such as the high-p-r azimuthal 
anisotropy [jj, the broadening and disappearance of di- 
jet correlations [ljj, extensions of the correlation analysis 
with respect to the reaction plane [ll| and possibilities for 
full jet and lost energy reconstruction |12| . The theory 
and phenomenology of multiparton dynamics in ultra- 
dense nuclear matter will, however, remain incomplete 
without a thorough investigation of the center of mass 
energy and system size dependence |l3l Il4j of medium 
induced non-Abelian gluon bremsstrahlung and the cor- 
responding hadron attenuation. Possible onset of pion 
suppression at the SPS energy of ^/snn = 17 GeV has 
been recently discussed . The intermediate y / ijv r Jv = 
62 GeV RHIC run is the next critical step in mapping 
out the elastic and inelastic scattering properties of a 
quark-gluon plasma via jet tomography [slllq. 

The magnitude of the energy loss driven nuclear 
quenching of moderate and high pt pion s is controlled 
by the soft parton rapidity density la.[lal. To relate the 
experimentally measured dN ch /dr] |ljfl to the effective 
dN 9 /dy we use \dr)/dy\ w 1.2 at r\ = y = 0, ignoring the 
^snn dependent changes in particle composition. The 
estimated effective 
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follows from the isospin symmetry of strong interac- 
tions and the approximate parton- hadron duality |18| . 



Straightforward application of Eq. (Q for central, 
Np ar t = 340, Au+Au collisions and dN ch /di] constrained 
from the data [l^ yields dN 9 /dy « 550, 850, 1150 at 
SPS, the intermediate and maximum RHIC energies, re- 
spectively. Such rapidity densities will likely be compat- 
ible with a soft participant scaling phenomenology 
The predicted snn — 200 GeV pion quenching 
was found to be in good agreement with the experimen- 
tal measurements p] but the original SPS tt° data j2fj 
would tend to disfavor parton energy loss in variance with 
the theoretical expectations Q . A more recent analysis of 
the low energy p + p baseline cross section shows that 
the WA98 [2JJ and CERES fH) data are not inconsis- 
tent with jet quenching calculations with dN 9 /dy = 400 
when the Cronin effect is taken into account 13] . The 
extracted effective gluon rapidity density however, still 
falls short of the expectation from the measured hadron 
multiplicities 0|. This deviation can be related to the 
uncertainties in the perturbative calculation in a theory 
with strong coupling and a non-negligible quark contribu- 
tion |23 to the bulk soft partons at ^/snn = 17 GeV at 
midrapidity. Therefore, the importance of studying the 
sensitivity of the observable spectral modification to a 
range of parton densities at a fixed center of mass energy 
should not be underestimated in theoretical calculations. 
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The purpose of this letter is to investigate the inter- 
play of the initial state multiple Cronin scattering and 
the final state energy loss in a thermalized QCD media 
with effective dN 9 /dy = 650 — 800 in moderate and large 
transverse momentum hadron production. Section II 
outlines the method for evaluating the medium induced 
gluon bremsstrahlung off fast partons in dense, dynam- 
ically expanding and finite nuclear matter. Section III 
presents the calculated depletion of pion multiplicities in 
central Au + Au reactions at y / sjv r /v = 62 GeV relative to 
the binary collision scaled p + p baseline. Summary and 
discussion is given in Section IV. 
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FIG. 1: (a) Medium induced single inclusive gluon spectrum dN g /dx and fractional intensity xdN g /dx versus x ~ co/E for a 
E = 10 GeV gluon jet. (b) Mean medium-induced gluon number (N g ) per jet for initial effective soft parton rapidity density 
dN 9 /dy = 650 — 800. (c) Corresponding mean fractional energy loss {AE)/E. (d) Probability distribution P(e) of the fractional 
energy loss of E = 10 GeV quarks and gluons via multiple independent gluon bremsstrahlung, e = 'Yl i '-^i/E. 



II. MEDIUM INDUCED GLUON 
BREMSSTRAHLUNG IN FINITE DYNAMICAL 
PLASMAS 

The full solution for the medium induced gluon ra- 
diation off jets produced in a hard collisions at early 
times Tj et ~ 1/E inside a nuclear medium of length L 



can be obtained to all orders in the correlations between 
the multiple scattering centers via the reaction operator 
approach 23] . Other existing techniques have been re- 
viewed in The double differential bremsstrahlung in- 
tensity for gluons with momentum k = [xp + , \i 2 /xp + , k] 
resulting from the sequential interactions of a fast parton 
with momentum p = [p + , 0, 0] can be written as 
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where ^2 = is understood. In the small angle eikonal propagators are denoted by 
limit x = k + /p + as u>/E. In Eq. J5J the color current 



C (m ,...,„) = iv k ln(k-q m q„) 2 

^(tn+l,— ,n)(,m,— ,n) ~ G(m+l ,n) ~ C( m ... >n ) . (3) 
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The momentum transfers are distributed according to 
a normalized elastic differential cross section, 



cr e l 
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which models scattering by soft partons with a ther- 
mally generated Debye screening mass fi(i). It has 
been shown [23L |24| via the cancellation of direct and 
virtual diagrams that in the eikonal limit only the 
gluon mean free path X g (i) enters the medium-induced 
bremsstrahlung spectrum in Eq. @. For gluon domi- 
nated bulk soft matter a e i(i) w §7ra 2 / p, 2 (i) and X g (i) = 
1 / a e i{i) p{i) . The characteristic path length dependence 
of the non-Abelian energy loss in Eq. @ comes from the 
interference phases and is differentially controlled by the 
inverse formation times, 



(k - q m q w ) 2 

2xE 
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and the separations of the subsequent scattering centers 
Azk = Zk — Zk—i- It is the non-Abelian analogue of the 
Landau-Pomeranchuk-Migdal destructive interference ef- 
fect in QED HI. 

For the case of local thermal equilibrium we relate 
all dimensional scales in the problem to the tempera- 
ture of the medium T(i), for example, p 2 (i) = 47ra s T(i) 2 
and the elastic scattering cross section indicated above. 
The gluon radiative spectrum is evaluated numerically 
in an ideal 1+1D Bjorken expanding plasma |26(, p(i) = 

Po(To/n) a , ei = £o(Vo/ T i) Q+1 ' 3 witn ol = 1 and v s be- 
ing the speed of sound. The scaling with proper time 
of the temperature, the Debye screening mass p(i) and 
the gluon mean free path X g (i) naturally follow. To lead- 
ing power, the initial equilibration time To cancels in the 
evaluation of the bremsstrahlung integrals Q since 



dn p(i)n 



1 dN 9 
A±_ dy 



(6) 



In Eq. © A± is the transverse size of the medium and 
dN 9 /dy is the relevant physical quantity that controls 
the attenuation of the final state partonic flux. 

Transverse, 0t 7^ 0, expansion was not explicitly in- 
cluded in the energy loss calculation. Hence, we first 
clarify its impact on the medium induced non-Abelian 
bremsstrahlung. Numerically, the contributions to the 
radiative spectrum in Eq. @ have to be evaluated in the 
background of the dynamically evolving soft parton mul- 
tiplicity. In practice, for realistic 3+1D hydrodynamic 
simulation only the dominant n — 1 term |23j has been 
considered in the mean energy loss approximation |27j . 
Analytic treatments of transverse expansion must there- 
fore provide important guidance to its effect on the ra- 
diative spectra and (AE) . 

The simplest approach to (3t 7^ would be to modify 
the power a of the 1+1D Bjorken case to emulate 3+ ID 
dynamics. Indeed, beyond leading power, Eq. iJBJ, a 2 (p) 



effectively lowers the value of Bjorken a relative to the 
naive fixed coupling result. It has also been argued [2|J 
that if the transport properties of the medium are re- 
lated to the energy density, deviations from the ideal 
plasma limit would lead to a < 1 via v 2 = (3 + A) -1 

with A = ^(^r) 2 + ■■■ • Technically, if a > 1 the 
density integrals that control the energy loss can still be 
defined 28] . However, for a medium of fixed size L where 
simple analytic results can be obtained, Eq. 10, or even 
for a spatially varying density profile po = /Oo( x t), such 
physical picture corresponds to a superluminal Bjorken 
expansion. In this scenario, the dilution of the quark- 
gluon plasma from the transverse flow of the soft partons 
is modeled by streaming along the collision axis. As a 
result, the partons are lost as potential scatterers for the 
hard jets that escape the plasma at y = and the energy 
loss is reduced. Another consequence of choosing a > 1 
is the decrease of the accumulated transverse momentum 
(fcy) oc J po(To/T) a dT with the size of the medium L 
or, in realistic A + A collisions, significantly lowering the 
dependence of (fef,) on centrality. This contradicts the 
measured steady growth of the di-jet acoplanarity versus 
Np ar t [T(i| which even exceeds the theoretical expecta- 
tions for the Bjorken a = 1 case ||. 

Varying a and attributing the deviation from unity, 
a — 1, to a model of spatially non-uniform medium leads 
to power law behavior p(xy) = Pq\xt — x^ol" 1 rather 
than the smooth Woods-Saxon dependence of the nuclear 
matter density. Additionally, there is a large ambiguity in 
the choice of a since as a function of the azimuthal angle 
4> of the jet propagation relative to the reaction plane 
the spatial density can either increase or decrease. We 
thus conclude that a > 1 is not a good emulation of the 
realistic (3t 7^ expansion and that the geometry profile 
and the soft parton dynamics should not be substituted 
for each other. 

An analytically tractable approximation that illus- 
trates the effects of transverse expansion has been de- 
veloped in [29} . In this scenario, for a medium of mean 
density po and mean radius R = L 



P(r) 



1 dN 9 1 
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The additional dilution relative to the 1+1D Bjorken case 
now generated by the transverse motion of the soft par- 
tons and the increasing transverse size A± of the system. 
For a discussion of non-central collisions, R x 7^ R y and 
p Tx ^ p Tyj S ee [23. It follows from Eq. Q that the 
propagating jets interact with the bulk matter at midra- 
pidity over an increased time period At ~ L/(l — (3t) > 
L that largely compensates for the rarefaction of the 
medium. It has been demonstrated [2^| that for small 
and moderate expansion velocities (3t the 1+3D angular 
((f)) averaged energy loss approximates well the 1+1D re- 
sult. Logarithmic corrections arise only in the (3r — > 1 
limit. It should now be physically intuitive why the 
Bjorken expansion is a good approximation in the cal- 
culation of quenching in the single and double inclusive 
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hadron spectra. 

An important aspect of the application of the theory 
of medium induced non-Abelian bremsstrahlung is the 
inclusion of finite kinematic bounds. These are partic- 
ularly relevant at RHIC and SPS energies over the full 
accessible px range. Additional details on the evaluation 
of the radiative energy loss are given in |23| . 

Figure 1(a) shows the differential gluon spectrum 
dNg/dx and the fractional radiation intensity xdN g /dx 
computed to 3-rd order in opacity (n = 1,2,3) from 
Eq. @ for a 10 GeV gluon jet. The fluctuations in the 
calculation reflect the numerical accuracy of the cancel- 
lation between the propagator poles, Eq. © , and the in- 
terference phases, Eq. (jSJ, and come predominantly from 
higher orders in opacity. These, however, do not affect 
the evaluation of the nuclear modification factor. The jet 
quenching strength is largely set by the n = 1 term [23j. 

In a thermalized medium the plasmon frequency 
0Jpi{i) ~ regulates the infrared modes and the 

medium induced bremsstrahlung results in a few semi- 
hard gluons. Numerical results for quark jets are not 
given since they differ by a simple Cf/Ca = 4/9 color 
factor. Figures 1(b) and 1(c) show the mean induced 
gluon number (N g ) and the mean fractional energy loss 
(AE)/E calculated directly from Eq. for a range of 
soft parton rapidity densities dN 9 /dy = 650 — 800. Com- 
paring the shape of (AE)/E in Figure 1(c) to the naive 
analytic expectation in the infinite kinematic limit 



(AE) 9C R ira 3 s 1 



E 



A, 



dN g 1 2E 
~d^ L E ln T^L 



-7T7 + --- (8) 



and observing the deviation from the hyperbolic l/E de- 
pendence it is easy to recognize the need for a careful 
treatment of phase space. It is also instructive to note 
that the average energy per gluon (AE) / (N g ) ~ 0.8 — 
1.9 GeV, which implies that the radiative quanta might 
be experimentally observable with finite px cuts 0, El • 
Applications that extend beyond the mean energy loss 
and invoke a probabilistic treatment with multiple gluon 
fluctuations require additional assumptions |30|. So far 
even the case of two gluon emission with scattering has 
not been calculated. The Poisson approximation assumes 
independent radiation and equivalence of the single inclu- 
sive and the single exclusive gluon spectrum. Finite kine- 
matics alone is guaranteed to violate this ansatz and so 
is angular ordering. Nevertheless, the usefulness of such 
probabilistic treatment is in allowing the system to max- 



imize the observable cross section and correspondingly 
minimize, to the extent to which this is possible, the ef- 
fect of the non-Abelian medium-induced bremsstrahlung. 
The probability P(e) for fractional energy loss e = 
LOi/E due to multiple gluon emission is given in Fig- 
ure 1(d). The e = bin contribution S(e)e~^ N <i' > is 
not shown. Details on the calculation of P(e) are given 
in |3(i| , where a fast iterative procedure for its evaluation 
form the gluon radiative spectrum dN g /dx was devel- 
oped. If a non-zero fraction of the probability P(e) falls 
in the e > 1 region, it is uniformly redistributed in the 
physical e € [0, 1] interval. This approach ensures that 
(AE) / E < 1 in the large energy loss limit and provides 
corrections relative to its direct evaluation from Eq. (0) 
that are different for quarks and gluons. The properties 
of P(e) can be summarized by the normalization of the 
first two e k , k = 0, 1 moments: 



/ P(e)de=l, [ eP(e)de 
Jo Jo 



E 



(9) 



We note the distinct difference in the manifestation of 
the large energy loss relative to the single inclusive spec- 
tra. In the former case the overall scale changes and in 
the probabilistic interpretation the distribution is shifted 
toward large e values, see Figure 1(d). 



III. PHENOMENOLOGICAL APPLICATION AT 
INTERMEDIATE ENERGIES 

Dynamical nuclear effects in p + A and A + A reactions 
are most readily detectable through the nuclear modifi- 
cation ratio 

h dN* A {h) T AA {h)da™ 

R ^ iPTl) - d^ll d yi d 2 p T1 mA + A > 

(10) 

where T AA (b) = / d 2 vT A (v)T B {r - b) is calculable 
in terms of the nuclear thickness functions T A {r) — 
J dz pa{y, z). In R^ a (pti) the uncertainty associated 
with the next-to leading order factor, Knlo, drops out. 

The standard lowest order perturbative expression for 
the hadron multiplicity in A + A reactions including the 
effects of multiple initial state scattering and the final 
state energy loss reads 0, : 
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FIG. 2: (a) Calculated nuclear modification factor Raa(pti) versus the center of mass energy y^ijvjv = 17, 62, 200 GeV for 
central Au + Au collisions. The yellow band represents the perturbative QCD expectation for the depletion of the neutral pion 
multiplicity at ^Jsnn = 62 GeV. Enhancement, arising from transverse momentum diffusion in cold nuclear matter without 
final state energy loss is given for comparison. PHENIX [|, WA98 [H and CERES [U data as presented in [ljj are shown, 
(b) The ratio of the attenuation of the single and double inclusive ty° cross sections Raa^Ptx) / Ra\ 2 (jP T i ~ Pt 2 ) m central 
Au+Au collisions at ^/snn = 62 GeV and (c) the manifestation of the double inclusive hadron suppression in the quenching 
of the away-side di-hadron correlation function C*2(A0) |39|l . 



In Eq. |(TT]| f a /A(x,Q 2 ) are the isospin corrected (A = 
Z + N) lowest order parton distribution functions [32| . 
S a / A {x, Q 2 ) is the leading twist shadowing parameteriza- 
tion |33| . and D(z)„/ r (z, Q 2 ) is the fragmentation func- 
tion into pions [34| . Vacuum and medium-induced ini- 
tial state parton broadening, (k^) = (ky) pp + (ky) nuc ;., 
is incorporated via a normalized Gaussian k^ smearing 
function [3l| . If the final state parton looses a fraction e 
of its energy the correspondingly rescaled fragmentation 
momentum fraction reads z* = 2/(1 — e). For consistency 
the calculation is performed in the same way as in [l3| 
where additional details can be found. 

In the pt < 5 GeV range, experimentally accessible 
at y / s]vAr = 62 GeV, sizable non-perturbative effects in 
baryon production, manifest in enhanced p/ir and A/if 
ratios, will likely be observed. Discussion of the mod- 
erate pt baryon phenomenology is bey ond the scope of 
this letter and details are given in [35j. In the limit of 
vanishing baryon masses, tub — > 0, perturbative calcula- 
tions of the nuclear modification factor ^^(pri) yield 
results comparable to the one for neutral and charged pi- 
ons. We finally note that the dominant contribution to 
the nuclear shadowin g may c ome from enhanced dynami- 
cal power corrections |36l |3jj , resulting from the multiple 
initial and final state interactions of the partons on a nu- 



cleus. However, in the calculated pr range (x > 0.8, 
Q 2 > 8 GeV 2 ) their effect was found to be small. 

Results form the perturbative calculation of the nu- 
clear modification to the neutral and charged pion pro- 
duction in central Au + Au collisions are given in the left 
hand side of Figure 2(a). At all energies there is a strong 
cancellation between the Cronin enhancement, which 
arises from the transverse momentum diffusion of fast 
partons in cold nuclear matter 0,H^|, and the subsequent 
inelastic jet attenuation in the final state. This interplay 
is most pronounced at the low SPS ^/snn = 17 GeV 
where, in the absence of quenching, the corresponding 
enhancement could reach a factor of 3-4 01 • With final 
state energy loss taken into account, R aa (pti) is shown 
versus the reanalyzed [HI WA98 [H and the CERES 
Pb + Pb and Pb + Au data. We note that even the 
previously extracted nuclear modification at the SPS al- 
lows for inelastic final state interactions in a medium of 
dNS/dy = 200 0. 

At present, there are strong indications that the 

highest nuclear matter density reached in the early stages 
(t = 0.6 fm) of the ^/sWn = 200 GeV central Au + Au 
collisions at RHIC may be on the order of 100 times cold 
nuclear matter density, e co id = 0.14 GeV/fm 3 . Perturba- 
tive calculations with a corresponding dN 9 /dy — 1150 



6 



are compatible with the measured 7r° attenuation 
as illustrated in Figure 2(a). Quenching ratios of sim- 
ilar magnitude, ^(pn): R'cpiPTi) ~ 0.2 — 0.25, have 
also been predicted and measured for inclusive charged 
hadrons at pr > 5 GeV |lj . It is interesting to note that 
in spite of the very large initial energy density the nu- 
clear attenuation is "only" a factor of 4 — 5. The reason 
for this somewhat unintuitive result is the strong longi- 
tudinal expansion in the absence of which energetic jets 
would have been completely absorbed. 

The yellow band represents a calculation of the ir° at- 
tenuation at y/sNN = 62 GeV from Eqs. (JTOJ, (JTTJ. The 
Cronin enhancement alone in central Au + Au reactions 
was found to be still sizable with ii^ (pri)max = 1-8 — 2 
at pt = 3 — 3.5 GeV and a subsequent decrease at 
higher transverse momenta. For comparison, at ^/snn — 
200 GeV the Cronin maximum decreases with the cen- 
ter of mass energy and R AA {pTi)max = 1.7. While such 
enhancement may naively appear large, it is consistent 
with the system (Au + Au versus d + Au) and ^s^n 
dependence coming from p T diffusion 0]. We note that 
even at the maximum RHIC energy the y = central 
d + Au collisions exhibit ~ 35% maximum pion enhance- 
ment |38| . The final state partonic energy loss was com- 
puted for a range of initial effective gluon rapidity den- 
sities dN 9 /dy = 650 — 800 as discussed in Section II. 
Figure 2(a) clearly shows that the perturbative calcula- 
tion presented here predicts a dominance of the inelastic 
jet interactions and net tt° quenching at the intermediate 
RHIC energy. 

Additional constraints for the energy loss calculations 
would arise from the measurement of a suppressed double 
inclusive hadron production 

Rhlh2 = dN AA (b) , T AA (b) do™ 
AA dyidy2d 2 pTid 2 pT2 dyidy2d 2 pTid 2 pT2 

(12) 

The cross sections in Eq. (|12f) can be evaluated as in |3|j 
by fragmenting the second hard scattered parton. We 
note that by unitarity the details of its energy loss and 
fragmentation do not affect the calculation of the sin- 
gle inclusive spectra. Di-hadron correlations, however, 
pick only one of the particle states. Since the second 
parent parton also looses a fraction of its energy in the 
medium, qualitatively 1 < R AA /R AA 2 < 2. Numer- 
ical estimates in Figure 2(b) from jet quenching alone 
indicate that the double inclusive hadron suppression is 
25% — 40% larger than the single inclusive quenching in 
the 5 < pti — Pti < 7 GeV range at ^snn = 62 GeV. 
The ratio was computed in the mean AE approximation, 
but with an explicit average over the di-jet production 
point for a realistic Woods-Saxon geometry. Double in- 
clusive cross section quenching is manifest in the atten- 
uation of the away-side correlation function C2(A0) = 
(l/N ttig )dN hxh2 /dA<f> [13. In Figure 2(c) it leads to a 
3 — 5 fold attenuation of the area A[ al relative to the 
p + p case. In contrast, elastic transverse momentum dif- 
fusion will only result in a broader C2(A0). We note 



that the experimentally measured suppression value will 
be sensitive to the subtraction of the elliptic flow V2 com- 
ponent |4l| , which may remove part or all of the residual 
correlations. The results shown in Figure 2(c), therefore, 
correspond to the smallest anticipated observable disap- 
pearance of the away-side jet at high px- 



IV. DISCUSSION AND SUMMARY 

For pr > 2 GeV the reduction of the pion multiplicity 
in central Au + Au reactions relative to the binary colli- 
sion scaled p + p result at the intermediate RHIC energy 
of ^snn = 62 GeV was found to be a factor of 2 — 3. 
This result is consistent with R aa (pt — 4 GeV) ^0.5 
from |40| . The moderate transverse momentum depen- 
dence of the quenching ratio arises from the interplay 
of the Cronin effect, in particular its decrease toward 
larger transverse momenta, and the jet energy loss. The 
medium induced energy loss also results in an apparent 
disappearance of the away-side jet at high px, which is 
similar to the ^/snn = 200 GeV result [13] • 

To relate the effective dN 9 jdy = 650 - 800 to the tem- 
perature and energy density we employ an initial equi- 
libration time To ~ 0.8 — 1 fm suggested by the hydro- 
dynamic description of the SPS data |4l|. For 1+1D 
expansion e = 6 - 8 GeV/fm 3 , T Q = 300 - 325 MeV and 
the lifetime of the plasma phase t — tq = 3.5 — 4 fm. Such 
initial conditions are already significantly above the cur- 
rent expectation for the critical temperature and energy 
density for a deconfinement phase transition, e c ~ 1 GeV 
and T c ~ 175 MeV ; 42j. 

From a theoretical point of view, the most interesting 
result would be a strong deviation of the pion attenu- 
ation from the pQCD calculation in Figure 2(a). This 
may be either a factor of ~ 5 suppression, comparable to 
the y/sNN = 200 GeV quenching, or binary and above 
binary scaling, as seen at the SPS energies. Both cases 
would indicate a marked non-linear dependence of the 
non-Abelian energy loss on the soft parton rapidity den- 
sity dN 9 jdy if we assume that parton- hadron duality still 
holds. Equally striking will be an absence of strong atten- 
uation in the away-side C2(A</>). It is important to check 
whether the established connection between the single 
and double inclusive hadron suppression persist at 
the intermediate RHIC energies, see Figure 2(c). 

In summary, the upcoming ^/Jnn = 62 GeV measure- 
ments and their comparison to theoretical calculations 
will provide a critical test of jet tomography. They will 
help to further clarify the properties of the hot and dense 
nuclear matter crated in relativistic nuclear collisions and 
the position that it occupies on the QCD phase transition 
diagram. 

Note added: Shortly after the completion of this work 
preliminary PHENIX data on ttq production and attenua- 
tion at the intermediate RHIC energy of ^snn — 62 GeV 
became available. Comparison between the data and the 
predictions given in this letter can be found in [43| . 
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